Apparent non-canonical trans-splicing is generated by reverse transcriptase in vitro

Trans-splicing, the in vivo joining of two RNA molecules, is well characterized in several groups of simple organisms but was long thought absent from fungi, plants and mammals. However, recent bioinformatic analyses of expressed sequence tag (EST) databases suggested widespread trans-splicing in mammals^1-2^. Splicing, including the characterised trans-splicing systems, involves conserved sequences at the splice junctions. Our analysis of a yeast non-coding RNA revealed that around 30% of the products of reverse transcription lacked an internal region of 117 nt, suggesting that the RNA was spliced. The junction sequences lacked canonical splice-sites but were flanked by direct repeats, and further analyses indicated that the apparent splicing actually arose because reverse transcriptase can switch templates during transcription^3^. Many newly identified, apparently trans-spliced, RNAs lacked canonical splice sites but were flanked by short regions of homology, leading us to question their authenticity. Here we report that all reported categories of non-canonical splicing could be replicated using an in vitro reverse transcription system with highly purified RNA substrates. We observed the reproducible occurrence of ostensible trans-splicing, exon shuffling and sense-antisense fusions. The latter generate apparent antisense non-coding RNAs, which are also reported to be abundant in humans^4^. Different reverse transcriptases can generate different products of template switching, providing a simple diagnostic. Many reported examples of splicing in the absence of canonical splicing signals may be artefacts of cDNA preparation

Association with Aurora-A controls N-MYC-dependent promoter escape and pause release of RNA polymerase II during the cell cycle

MYC proteins bind globally to active promoters and promote transcriptional elongation by RNA polymerase II (Pol II). To identify effector proteins that mediate this function, we performed mass spectrometry on N-MYC complexes in neuroblastoma cells. The analysis shows that N-MYC forms complexes with TFIIIC, TOP2A, and RAD21, a subunit of cohesin. N-MYC and TFIIIC bind to overlapping sites in thousands of Pol II promoters and intergenic regions. TFIIIC promotes association of RAD21 with N-MYC target sites and is required for N-MYC-dependent promoter escape and pause release of Pol II. Aurora-A competes with binding of TFIIIC and RAD21 to N-MYC in vitro and antagonizes association of TOP2A, TFIIIC, and RAD21 with N-MYC during S phase, blocking N-MYC-dependent release of Pol II from the promoter. Inhibition of Aurora-A in S phase restores RAD21 and TFIIIC binding to chromatin and partially restores N-MYC-dependent transcriptional elongation. We propose that complex formation with Aurora-A controls N-MYC function during the cell cycle

Characterization and simulation of cDNA microarray spots using a novel mathematical model

<p>Abstract</p> <p>Background</p> <p>The quality of cDNA microarray data is crucial for expanding its application to other research areas, such as the study of gene regulatory networks. Despite the fact that a number of algorithms have been suggested to increase the accuracy of microarray gene expression data, it is necessary to obtain reliable microarray images by improving wet-lab experiments. As the first step of a cDNA microarray experiment, spotting cDNA probes is critical to determining the quality of spot images.</p> <p>Results</p> <p>We developed a governing equation of cDNA deposition during evaporation of a drop in the microarray spotting process. The governing equation included four parameters: the surface site density on the support, the extrapolated equilibrium constant for the binding of cDNA molecules with surface sites on glass slides, the macromolecular interaction factor, and the volume constant of a drop of cDNA solution. We simulated cDNA deposition from the single model equation by varying the value of the parameters. The morphology of the resulting cDNA deposit can be classified into three types: a doughnut shape, a peak shape, and a volcano shape. The spot morphology can be changed into a flat shape by varying the experimental conditions while considering the parameters of the governing equation of cDNA deposition. The four parameters were estimated by fitting the governing equation to the real microarray images. With the results of the simulation and the parameter estimation, the phenomenon of the formation of cDNA deposits in each type was investigated.</p> <p>Conclusion</p> <p>This study explains how various spot shapes can exist and suggests which parameters are to be adjusted for obtaining a good spot. This system is able to explore the cDNA microarray spotting process in a predictable, manageable and descriptive manner. We hope it can provide a way to predict the incidents that can occur during a real cDNA microarray experiment, and produce useful data for several research applications involving cDNA microarrays.</p

Gene expression profiling of gliomas: merging genomic and histopathological classification for personalised therapy

The development of DNA microarray technologies over the past decade has revolutionised translational cancer research. These technologies were originally hailed as more objective, comprehensive replacements for traditional histopathological cancer classification systems, based on microscopic morphology. Although DNA microarray-based gene expression profiling (GEP) remains unlikely in the near term to completely replace morphological classification of primary brain tumours, specifically the diffuse gliomas, GEP has confirmed that significant molecular heterogeneity exists within the various morphologically defined gliomas, particularly glioblastoma (GBM). Herein, we provide a 10-year progress report on human glioma GEP, with focus on development of clinical diagnostic tests to identify molecular subtypes, uniquely responsive to adjuvant therapies. Such progress may lead to a more precise classification system that accurately reflects the cellular, genetic, and molecular basis of gliomagenesis, a prerequisite for identifying subsets uniquely responsive to specific adjuvant therapies, and ultimately in achieving individualised clinical care of glioma patients

Radiogenomic Mapping of Edema/Cellular Invasion MRI-Phenotypes in Glioblastoma Multiforme

Despite recent discoveries of new molecular targets and pathways, the search for an effective therapy for Glioblastoma Multiforme (GBM) continues. A newly emerged field, radiogenomics, links gene expression profiles with MRI phenotypes. MRI-FLAIR is a noninvasive diagnostic modality and was previously found to correlate with cellular invasion in GBM. Thus, our radiogenomic screen has the potential to reveal novel molecular determinants of invasion. Here, we present the first comprehensive radiogenomic analysis using quantitative MRI volumetrics and large-scale gene- and microRNA expression profiling in GBM.Based on The Cancer Genome Atlas (TCGA), discovery and validation sets with gene, microRNA, and quantitative MR-imaging data were created. Top concordant genes and microRNAs correlated with high FLAIR volumes from both sets were further characterized by Kaplan Meier survival statistics, microRNA-gene correlation analyses, and GBM molecular subtype-specific distribution.The top upregulated gene in both the discovery (4 fold) and validation (11 fold) sets was PERIOSTIN (POSTN). The top downregulated microRNA in both sets was miR-219, which is predicted to bind to POSTN. Kaplan Meier analysis demonstrated that above median expression of POSTN resulted in significantly decreased survival and shorter time to disease progression (P<0.001). High POSTN and low miR-219 expression were significantly associated with the mesenchymal GBM subtype (P<0.0001).Here, we propose a novel diagnostic method to screen for molecular cancer subtypes and genomic correlates of cellular invasion. Our findings also have potential therapeutic significance since successful molecular inhibition of invasion will improve therapy and patient survival in GBM

Functional microarray analysis suggests repressed cell-cell signaling and cell survival-related modules inhibit progression of head and neck squamous cell carcinoma

<p>Abstract</p> <p>Background</p> <p>Cancer shows a great diversity in its clinical behavior which cannot be easily predicted using the currently available clinical or pathological markers. The identification of pathways associated with lymph node metastasis (N+) and recurrent head and neck squamous cell carcinoma (HNSCC) may increase our understanding of the complex biology of this disease.</p> <p>Methods</p> <p>Tumor samples were obtained from untreated HNSCC patients undergoing surgery. Patients were classified according to pathologic lymph node status (positive or negative) or tumor recurrence (recurrent or non-recurrent tumor) after treatment (surgery with neck dissection followed by radiotherapy). Using microarray gene expression, we screened tumor samples according to modules comprised by genes in the same pathway or functional category.</p> <p>Results</p> <p>The most frequent alterations were the repression of modules in negative lymph node (N0) and in non-recurrent tumors rather than induction of modules in N+ or in recurrent tumors. N0 tumors showed repression of modules that contain cell survival genes and in non-recurrent tumors cell-cell signaling and extracellular region modules were repressed.</p> <p>Conclusions</p> <p>The repression of modules that contain cell survival genes in N0 tumors reinforces the important role that apoptosis plays in the regulation of metastasis. In addition, because tumor samples used here were not microdissected, tumor gene expression data are represented together with the stroma, which may reveal signaling between the microenvironment and tumor cells. For instance, in non-recurrent tumors, extracellular region module was repressed, indicating that the stroma and tumor cells may have fewer interactions, which disable metastasis development. Finally, the genes highlighted in our analysis can be implicated in more than one pathway or characteristic, suggesting that therapeutic approaches to prevent tumor progression should target more than one gene or pathway, specially apoptosis and interactions between tumor cells and the stroma.</p

Large-scale integration of cancer microarray data identifies a robust common cancer signature

<p>Abstract</p> <p>Background</p> <p>There is a continuing need to develop molecular diagnostic tools which complement histopathologic examination to increase the accuracy of cancer diagnosis. DNA microarrays provide a means for measuring gene expression signatures which can then be used as components of genomic-based diagnostic tests to determine the presence of cancer.</p> <p>Results</p> <p>In this study, we collect and integrate ~ 1500 microarray gene expression profiles from 26 published cancer data sets across 21 major human cancer types. We then apply a statistical method, referred to as the <it>T</it>op-<it>S</it>coring <it>P</it>air of <it>G</it>roups (TSPG) classifier, and a repeated random sampling strategy to the integrated training data sets and identify a common cancer signature consisting of 46 genes. These 46 genes are naturally divided into two distinct groups; those in one group are typically expressed less than those in the other group for cancer tissues. Given a new expression profile, the classifier discriminates cancer from normal tissues by ranking the expression values of the 46 genes in the cancer signature and comparing the average ranks of the two groups. This signature is then validated by applying this decision rule to independent test data.</p> <p>Conclusion</p> <p>By combining the TSPG method and repeated random sampling, a robust common cancer signature has been identified from large-scale microarray data integration. Upon further validation, this signature may be useful as a robust and objective diagnostic test for cancer.</p

Deep RNA sequencing analysis of readthrough gene fusions in human prostate adenocarcinoma and reference samples

<p>Abstract</p> <p>Background</p> <p>Readthrough fusions across adjacent genes in the genome, or transcription-induced chimeras (TICs), have been estimated using expressed sequence tag (EST) libraries to involve 4-6% of all genes. Deep transcriptional sequencing (RNA-Seq) now makes it possible to study the occurrence and expression levels of TICs in individual samples across the genome.</p> <p>Methods</p> <p>We performed single-end RNA-Seq on three human prostate adenocarcinoma samples and their corresponding normal tissues, as well as brain and universal reference samples. We developed two bioinformatics methods to specifically identify TIC events: a targeted alignment method using artificial exon-exon junctions within 200,000 bp from adjacent genes, and genomic alignment allowing splicing within individual reads. We performed further experimental verification and characterization of selected TIC and fusion events using quantitative RT-PCR and comparative genomic hybridization microarrays.</p> <p>Results</p> <p>Targeted alignment against artificial exon-exon junctions yielded 339 distinct TIC events, including 32 gene pairs with multiple isoforms. The false discovery rate was estimated to be 1.5%. Spliced alignment to the genome was less sensitive, finding only 18% of those found by targeted alignment in 33-nt reads and 59% of those in 50-nt reads. However, spliced alignment revealed 30 cases of TICs with intervening exons, in addition to distant inversions, scrambled genes, and translocations. Our findings increase the catalog of observed TIC gene pairs by 66%.</p> <p>We verified 6 of 6 predicted TICs in all prostate samples, and 2 of 5 predicted novel distant gene fusions, both private events among 54 prostate tumor samples tested. Expression of TICs correlates with that of the upstream gene, which can explain the prostate-specific pattern of some TIC events and the restriction of the <it>SLC45A3-ELK4 </it>e4-e2 TIC to <it>ERG</it>-negative prostate samples, as confirmed in 20 matched prostate tumor and normal samples and 9 lung cancer cell lines.</p> <p>Conclusions</p> <p>Deep transcriptional sequencing and analysis with targeted and spliced alignment methods can effectively identify TIC events across the genome in individual tissues. Prostate and reference samples exhibit a wide range of TIC events, involving more genes than estimated previously using ESTs. Tissue specificity of TIC events is correlated with expression patterns of the upstream gene. Some TIC events, such as <it>MSMB-NCOA4</it>, may play functional roles in cancer.</p